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Abstract 


A  rapid  ship  hull  form  generation  tool  is  developed  in  support  of  an  ongoing  ONR 
research  project  on  the  development  of  a  methodology  for  innovative  hydrodynamic 
design  of  ship  hull  forms.  The  hull  lines/curves  in  this  tool  are  described  by  various 
orders  of polynomials  and  the  polynomials  are  determined  using  hull  form  parameters. 
The  hull  generation  tool  is  implemented  using  Matlab  to  enable  graphic  user  interface. 
Two  hull  forms  are  used  for  the  validation.  The  present  tool  is  able  to  generate  the  same 
hull  forms  as  the  sample  hull  forms  given  by  Bjorklund  &  Fuller  (1976)  with  given 
parameters.  This  rapid  hull  form  generation  tool  describes  hull  form  using  naval 
architect ’s  language;  generates  hull  form  from  scratch  in  terms  of  hull  form  parameters; 
and  establishes  the  link  between  the  hull  form  and  form  parameters.  It  can  be  used  to 
generate  the  hull  form  for  hydrostatic  analysis  and  hydrodynamic  design  optimization  in 
the  early  stage  of  ship  design. 
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Introduction 

Objective 

In  order  to  carry  out  hydrodynamic  performance  analysis  and  design  optimization  in  the 
early  stage  of  ship  design,  it  is  essential  to  have  a  rapid  ship  hull  form  generation  tool  that 
can  be  used  to  generate  hull  forms  in  terms  of  design  needs  and  hull  form  parameters. 

The  objective  of  this  project  is  to  develop  a  rapid  hull  form  generation  tool  that  describes 
hull  form  using  naval  architect’s  language;  generates  a  hull  form  from  scratch  in  terms  of 
hull  form  parameters;  and  establishes  the  link  between  the  hull  form  and  form 
parameters.  The  tool  can  be  used  to  define  the  hull  form  for  hydrostatic  analysis  and 
hydrodynamic  design  optimization  in  the  early  stage  of  ship  design. 

Background 

Hydrodynamic  design  of  ships  involves  several  stages,  from  early-stage  and  preliminary 
design  to  final  design.  It  has  become  increasingly  important  to  evaluate  multiple 
hydrodynamic  aspects  of  performance  simultaneously  and  efficiently  during  the  early 
stage  of  the  design  process. 

The  Computational  Fluid  Dynamics  (CFD)  group  at  George  Mason  University  has  an 
ongoing  research  project  sponsored  by  the  Office  of  Naval  Research  (ONR)  on  the 
development  of  a  methodology  for  innovative  hydrodynamic  design  of  ship  hull  forms. 
The  methodology  under  development  will  lead  to  a  design  tool  to  explore  both 
conventional  and  non-traditional  hull  forms,  to  define  hull  forms  that  exhibit  both 
superior  seakeeping  and  low  resistance  at  early  stage  of  ship  design,  and  to  investigate 
innovative  new  hull  forms  with  specific  design  constraints. 

In  order  to  apply  CFD-based  tool  to  the  hydrodynamic  design  optimization  of  ship  hull 
forms,  an  initial  hull  form  is  required.  The  hull  surface  can  then  be  modified  in  terms  of 
the  design  parameters  associated  with  the  hull  surface  representation  and  modification 
techniques.  Two  hull  form  modification  techniques  have  been  developed  to  obtain 
maximal  benefits  from  both  conventional  and  parametric  modeling  techniques  (Yang  et 
al.,  2008;  Kim  et  al.  2010): 

(i)  A  modified  NURBS  technique  is  combined  with  a  parametric  global  hull 
modification  technique  by  varying  the  sectional  area  curve  with  a  shifting  method, 
in  which  the  design  variables  in  the  NURBS  technique  can  be  reduced  via  a 
grouping  method  after  introducing  dependent  and  independent  design  variables; 

(ii)  A  radial  basis  function  (RBF)  approach  is  combined  with  a  parametric  global  hull 
modification  technique  by  varying  the  sectional  area  curve  with  a  shifting  method. 
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The  first  technique  allows  for  larger  variation,  and  the  second  one  permits  the  constraints, 
such  as  fixed  waterline  or  section  shape,  to  be  implemented  easily.  Figs.  1.1  and  1.2  show 
the  preliminary  optimization  results  of  the  Series  60  hull  model  by  minimizing  the  wave 
resistance  using  technique  (i)  for  the  hull  form  modification  in  the  optimization  process. 
As  it  is  shown  in  Fig.  1.1,  a  new  stem  profile  is  generated.  Figs.  2.1  and  2.2  show  the 
preliminary  bow-region  optimization  results  of  the  Model  5415  hull  by  minimizing  the 
wave  resistance  using  technique  (ii)  for  the  hull  form  modification  in  the  optimization 
process.  As  it  is  shown  in  Fig.  2.1,  new  sections  and  profiles  are  generated.  It  can  be 
observed  from  Fig.  1.2  and  Fig.  2.2  that  both  optimal  hulls  can  yield  appreciable  wave 
drag  reduction. 


Initial  Hull:  Series  60  Hull  Model 


Fig.  1.1  Comparison  of  the  ship  lines  Fig.  1.2  Wave  drag  reductions  at 

(sections  &  profiles)  the  design  speed  range 


Initial  Hull:  DTMB  Model  5415  Hull 
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Fig.  2. 1  Comparison  of  the  ship  lines  (sections  Fig.  2.2  Wave  drag  reductions  at  the  design 
&  profiles)  speed  range 

These  two  combined  techniques  for  hull  form  modifications  have  been  proven  to  be  very 
efficient  and  effective.  However,  these  two  techniques  would  work  only  in  the  case  that 
an  initial  hull  form  exists.  It  is  of  great  interest  for  naval  architects  to  have  a  rapid  hull 
form  generation  tool  to  generate  the  hull  form  from  scratch  using  their  language,  namely 
in  terms  of  form  parameters,  and  link  this  tool  with  the  CFD-based  hull  form 
optimization  tool.  An  approach  given  by  Bjorklund  &  Fuller  (1976)  is  adopted  in  this 
project  to  describe  the  hull  lines/curves  by  various  orders  of  polynomials.  The 
polynomials  are  determined  using  hull  form  parameters.  The  following  sections  will 
discuss  the  formulation  and  implementation  of  the  hull  generation  tool  (HullGen_201 1) 
developed  in  this  project. 


HullGen  2011 


A  rapid  ship  hull  form  generation  tool,  HullGen_201 1  is  implemented  using  Matlab  to 
enable  the  graphic  user  interface.  An  approach  given  by  Bjorklund  &  Fuller  (1976)  is 
adopted  to  generate  the  hull  form  in  terms  of  the  hull  parameters.  The  hull  lines/curves 
are  described  by  polynomials  of  various  orders.  The  unknown  coefficients  of  the 
polynomials  are  determined  by  the  input  hull  form  parameters. 

The  hull  form  (body  plan)  in  HullGen  is  defined  in  terms  of  the  following  control  curves: 

•  Section  area  curve 

•  Load  waterline  curve 

•  Deck  edge  curve 

•  Sheer  profile 

•  Section  area  coefficient  curve 
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•  Keel  rise  curve 

•  Section  slope  curves 

•  Flat  bottom  curve 

Each  curve  is  defined  by  a  polynomial  with  a  given  order.  The  polynomials  are  described 
in  a  non-dimensional  coordinate  system,  which  implies  that  the  coordinates  of  the 
polynomials  are  normalized  using  different  hull  parameters.  Specifically,  x-coordinates 
(the  longitudinal  coordinates)  are  normalized  using  14  LPP,  the  length  between  forward 
and  aft  perpendiculars.  The  y-coordinates  are  normalized  using  different  parameters, 
depending  on  which  curve  is  being  considered. 

The  computer  code  HullGen_201 1  is  written  in  Matlab.  It  consists  of  a  main  script  and  a 
number  of  functions.  The  user  has  the  option  of  either  using  the  default  hull  form 
parameters  given  in  the  code,  or  defining  them  in  an  input  data  file.  The  user  can  then 
modify  the  hull  parameters  interactively  to  generate  various  curves/lines  and  obtain  hull 
body  plan.  Two  sample  input  data  files  are  listed  in  the  appendix,  which  can  be  used  to 
generate  an  aircraft  carrier  type  hull  form  and  a  destroyer/cruiser  type  hull  form, 
respectively.  The  user  can  type  number  “1”  or  “2”  in  the  following  prompt  box  to  input 
the  data  file  or  use  default  hull  form  parameters  after  the  code  HullGen  starts  to  run: 


Press  1  to  input  a  file  or  2  to  use  default  parameters? 


I  QK  J  ;CancelJ 


A  new  prompt  box  will  then  be  displayed  for  the  user  to  specify  the  name  of  the  input 
data  file  if  the  option  2  is  selected. 

The  mathematical  formulation  of  the  control  curves/hull  lines  and  the  user  graphic 
interface  prompt  boxes  displayed  with  the  default  hull  form  parameters  are  discussed  in 
the  following  sections. 

Section  Area  Curve 


After  the  user  decides  to  use  either  default  hull  form  parameters  given  by  the  computer 
code  or  the  hull  form  parameters  specified  in  the  input  data  file,  HullGen  will  ask  for  the 
parameters  shown  in  the  following  prompt  box  to  generate  a  section  area  curve: 
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^  ^  Section  Area  Curve 


Area  at  Station  0  (AO) 

0.001 

S  lope  of  Curve  at  Station  0  (S  O) 

-0.5 

Location  of  Desired  Station  of  Maximum  Area  (Xmax) 
-0.1 

S  lope  of  Curve  at  Xmax 
0 


Length  of  Parallel  Sections  (Lmid) 
0.001 

Area  at  Station  20 

0.05 


of  Curve  at  Station  20 

Prismatic  Coeff.  (Cp) 

063 

Longitudinal  Center  of  Buoyancy  (LCB) 

^3 

|  OK  J  |  Cancel  j 


Fig.  3.  Prompt  box  for  section  area  curve 

Hereafter,  all  the  graphic  user  interface  prompt  box  displays  the  default  hull  form 
parameters  given  in  the  code. 

There  are  two  options  to  consider  for  describing  the  section  area  curve,  one  of  which  is 
for  the  hull  form  without  a  parallel  mid  body,  the  other  with  a  parallel  mid  body.  The 
coordinates  in  the  section  area  curve  are  non-dimensionalized  such  that  the  distance  from 
the  forward  perpendicular  (Station  0)  to  the  aft  perpendicular  (Station  20)  makes  x  run 
from  +1.0  to  -1.0  with  midships  (Station  10)  at  jc=0.  y= 1  is  the  ordinate  at  the  station  of 
maximum  area. 


The  following  boundary  conditions  are  needed  to  define  the  section  area  curve  for  the 
hull  form  without  a  parallel  mid  body: 


(1) 

x  =  1.0,  y  =  Area0 

(Area  at  Station  0) 

(2) 

x  =  1.0,  y'  =  Slope0 

(Slope  of  the  area  curve  at  Station  0) 

(3) 

*  =  *ma  x>T  =  1-° 

(Location  of  the  station  of  maximum  area) 

(4) 

x  =  xBUX,y'  =  0.0 

(Slope  at  the  station  of  maximum  area) 

(5) 

x  =  -1.0,  y  =  Area20 

(Area  at  Station  20) 

(6) 

x  =  -1.0,  y'  =  Slope20 

(Slope  of  the  area  curve  at  Station  20) 

(7) 

j!,  ydx  =  2.0  *  Cp 

(Area  under  the  curve  defined  by  the 
prismatic  coefficient  CP ) 

(8) 

ill  yxdx  =  2.0  *Cp*  LCB 

(Centroid  of  the  area  under  the  curve, 

defined  by  the  longitudinal  center  of 
buoyancy  LCB) 

A  7th 

order  polynomial,  written  as  follows,  is  used  to  produce  the  section  area  curve: 
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y  —  Aq  +  A^x  +  A2x  +  A2x^  +  A^x  +  A$x^  +  A^x^  +  A2x ^ 

in  which  the  unknown  coefficients  Aq  to  A  7  are  determined  by  the  boundary  conditions 
(1)  to  (8)  described  above. 


If  the  user  wants  to  produce  a  hull  with  a  parallel  mid  body,  a  slightly  different  set  of 
boundary  conditions  is  necessary: 


(1)  x  =  1 .0,  y  =  Area0 

(2)  x  =  1 .0,  y'  =  Slope0 

(3)  x  =  x2  =  xmax  +  Lmid  /2,  y  =  1.0 

(4)  x  =  x2  =  xmax  +  Lmid  / 2,  /  =  0.0 

(5)  x  =  xj  =  xmax  -  Lmid 12,  y  =  1 .0 

(6)  X  =  X!  =  Xmax  -  Lmid  /  2,  /  =  0.0 

(7)  x  =  —1 .0,  y  =  Area20 

(8)  x  =  -1.0 ,  y' =  Slope20 

(9)  }_"j  ydx  +  Lmid  *1.0  +  J*2  ydx  =  2.0  *CP 

( 1 0)  j:v;  yxdx  +  Lmid  *  1 .0  *  xmax  +  \l2  yxdx  = 


(Area  at  Station  0) 

(Slope  of  the  area  curve  at  Station  0) 
(Location  of  the  forward  end  of  the  flat 
section  in  the  area  curve) 

(Slope  of  the  flat  section  of  area  curve) 

(Location  of  the  after  end  of  the  flat  section 
in  the  area  curve) 

(Slope  of  the  flat  section  of  area  curve) 
(Area  at  Station  20) 

(Slope  of  the  area  curve  at  Station  20) 

(Area  under  the  curve,  defined  by  the 
prismatic  coefficient  CP) 

2.0  *  CP*  LCB 

(Centroid  of  the  area  under  the  curve, 
defined  by  the  longitudinal  center  of 
buoyancy  LCB) 


The  above  conditions  can  be  used  to  generate  a  9th  order  polynomial.  The  section  area 
curve  is  represented  by  the  9th  order  polynomial  except  the  parallel  mid  body  part  with  a 
horizontal  line  of  value  1.0.  Fig.  4  shows  the  hull  section  area  curve  with  a  parallel  mid 
body  generated  using  the  parameters  displayed  in  Fig.  3. 


Fig.  4.  Hull  section  area  curve  with  a  parallel  mid  body 
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Load  Waterline  Curve 


The  load  waterline  curve  is  a  representation  of  the  shape  of  the  hull  plan  form  at  the 
design  waterline.  The  creation  of  this  curve  is  almost  identical  to  that  of  the  section  area 
curve,  except  in  place  of  the  prismatic  coefficient  and  the  longitudinal  center  of 
buoyancy,  the  load  waterline  curve  uses  the  waterplane  coefficient  and  the  longitudinal 
center  of  floatation,  respectively.  The  y-coordinates  are  normalized  using  one  half  of  the 
ship  beam  defined  at  the  station  of  the  maximum  area. 

After  HullGen  generates  the  section  area  curve,  it  will  ask  for  the  parameters  shown  in 
the  following  prompt  box  to  generate  a  load  waterline  curve: 


Load  Waterline  Curve 


Offset  at  Station  0(Y0) 
O.GQ5 


or  Curve  at  Station  0  (SO) 


Location  or  Desired  Station  or  Maximum  Area 


,xm“  - 


F 


or  Curve  at  Xmax 


Length  or  Parallel  Sections  (Lmid) 


0.07 


Offset  at  Station  20 


0.4063 


S  lope  or  Curve  at  Station  20 
1.65 


Waterplane  Coeff.  (Cvvp) 

0.7514 

Longitudinal  Center  or  Floation  (LCF) 


-6.0 


(  OK  J  |  Cancel  I 


Fig.  5.  Prompt  box  for  load  waterline  curve 


As  before,  there  are  two  options  to  consider.  If  the  user  does  not  want  a  parallel  mid 
body,  the  following  are  the  boundary  conditions  that  the  curve  should  satisfy: 

(1)  x  =  1 .0,  y  =  Y0  (Offset  at  Station  0) 

(2)  x  =  1 .0,  y'  =  Slope0  (Slope  of  the  load  waterline  curve  at 

Station  0) 

(3)  x  =  xmax ,  y  =  1.0  (Offset  at  xmax  and  xmax  is  the  station  of 

maximum  area  from  the  section  area  curve. 
This  is  where  the  beam  of  the  ship  is 
defined.) 

(4)  x  =  xmax,  y'  =  Slope  _xmax  (Slope  at  xmax  which  can  be  defined  so  that 

the  beam  on  the  design  waterline  can  be 
maximum  forward  of  xmax  ,  at  xmax  ,  or  aft  of 
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xmax  •  The  result  is  that  the  actual  beam  of 
the  ship  can  be  greater  than  the  beam  at  the 
station  of  maximum  area) 


(5) 

x  =  -1.0  ,y  =  Y20 

(Offset  at  Station  20) 

(6) 

x  =  -1.0,  y  =  Slope20 

(Slope  of  the  load  waterline  curve  at 
Station  20) 

(V) 

k 

* 

o 

<N 

II 

■S 

(Area  under  the  curve,  defined  by  the 
waterplane  coefficient  CWp) 

(8) 

J*  j  yxdx  =  2.0  *  CWP  *  LCF 

(Centroid  of  the  area  under  the  curve, 

defined  by  longitudinal  center  of 
floatation  LCF) 


The  above  boundary  conditions  can  be  used  to  create  a  7th  order  polynomial  to  represent 
the  load  waterline  curve  for  the  hull  that  does  not  have  a  parallel  mid  body.  However,  if 
the  user  wants  to  design  a  hull  with  a  parallel  mid  body,  a  slightly  different  set  of 
boundary  conditions  is  required: 


(1) 

x  =  1.0  ,y  =  Y0 

(Offset  at  Station  0) 

(2) 

x  =  1.0  ,y'  =  Slope0 

(Slope  of  the  load  waterline  curve  at 
Station  0) 

(3) 

X  =  X2  =*max+4«rf/2,y  =  1.0 

(Offset  at  forward  end  of  parallel  section) 

(4) 

*  =  *2  =*max+4,irf/2,/  =  0.0 

(Slope  of  parallel  section) 

(5) 

x  =  xl  =*max-4«rf/2,y  =  1.0 

(Offset  at  after  end  of  parallel  section) 

(6) 

x  =  xl  =  *max  -Kid1 2,  /  =  0.0 

(Slope  of  parallel  section) 

(7) 

x  =  -1.0  ,y  =  Y20 

(Offset  at  Station  20) 

(8) 

x  =  -1.0,  y  =  Slope20 

(Slope  of  the  load  waterline  curve  at 
Station  20) 

(9) 

\-\ ydx  +  Kid  * 1  -0  +  j'2  ydx  =  2.0  *  CWF 

,  (Area  under  the  curve,  defined  by  the 
waterplane  coefficient  Cwp ) 

(10) 

1-1  yxdx  +  Lmid  *  1 .0  *  xmax  +  flX2  yxdx  = 

2.0  *  CWP  *  LCF 

(Centroid  of  area  under  the  curve,  defined 

by  the  longitudinal  center  of  floatation  LCF) 

The  above  boundary  conditions  can  be  used  to  generate  a  9th  order  polynomial.  The  load 
waterline  curve  is  represented  by  the  9th  order  polynomial  except  the  parallel  mid  body 
section  of  the  curve  with  a  horizontal  line  of  value  1 .0.  Fig.  6  shows  the  load  waterline 
curve  with  a  parallel  mid  body  generated  using  the  parameters  displayed  in  Fig.  5. 
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Fig.  6  Load  waterline  curve  with  a  parallel  mid  body 


Deck  at  Edge  Curve 


The  deck  at  edge  curve  represents  the  shape  of  the  deck.  This  curve  is  created  by 
combining  three  curves:  i)  a  3ld  order  polynomial  for  the  curve  forward  of  a  flat  middle 
section;  ii)  a  straight  line  for  the  middle  section  if  any;  and  iii)  a  3rd  order  polynomial  for 
the  curve  aft  of  the  flat  middle  section.  The  y-coordinates  are  normalized  using  half  of  the 
ship  beam  defined  at  the  station  of  the  maximum  area. 


After  the  HullGen  generates  the  load  waterline  curve,  it  will  ask  for  the  parameters 
shown  in  the  following  prompt  box  to  generate  a  flush  deck  at  edge  curve: 


Flush  Deck  at  Edge  Curve 

Offset  at  Station  0(Y0) 

0.75 

S  lope  of  Curve  at  Station  0  (S  0) 

-0.25 

Location  of  Flat  Middle  Section  (Xmid) 

-0.1025 

^^3 

Length  of  Flat  Middle  Sections  (Lmid) 

0.22 

Offset  at  Forward  End  of  Flat  Middle  Section  (Yrnid) 

1.0 

Slope  of  Curve  at  Forward  End  of  Flat  Middle  Section  (Smid) 

0 

Offset  at  Station  20fY20) 

0.0 

S  lope  of  Curve  at  Station  20  ( S  20) 

0.35 

’ 

OK 

]  [Cancel] 

Fig.7.  The  prompt  box  for  the  flush  deck  at  edge  curve 
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To  create  the  forward  section  of  the  flush  deck  at  edge  curve,  the  following  boundary 
conditions  must  be  satisfied: 


(1) 

x  =  1.0  ,y  =  Y0 

(Offset  at  Station  0) 

(2) 

X 

II 

o 

X 

II 

9s 

c? 

(Slope  of  the  flush  deck  at  edge  curve  at 
Station  0) 

(3) 

x  -  x2 ,  y  -  Ymid 

(Offset  at  forward  end  of  flat  middle 
section  defined  by  TMid  and  xmid) 

(4) 

x  =  x2,  y'  =  0 

(Zero  slope  by  definition) 

rr\ 

The  above  conditions  can  be  used  to  generate  a  3  order  polynomial  for  the  forward 
section.  Note  how  there  are  no  coefficients  to  govern  the  shape  of  the  curve.  The  aft 
portion  of  the  curve  is  formulated  in  a  similar  manner  as  above.  Fig.  8  shows  a  flush  deck 
at  edge  curve  generated  using  the  parameters  displayed  in  Fig.  7. 


Fig.  8.  Flush  deck  at  edge  curve 


Sheer  Profile 


The  sheer  profile  defines  the  depths  at  various  points,  which  also  defines  the  heights  of 
the  deck  at  every  section.  To  generate  this  curve,  the  user  inputs  the  depths  at  Stations  0, 
10,  and  20.  With  these  boundary  conditions,  HullGen  generates  a  2nd  order  polynomial. 
This  curve  is  dimensionalized,  i.e.  the  user  inputs  the  actual  values  of  the  depths  instead 
of  relative  values.  Fig.  9  shows  the  prompt  box  for  the  sheer  profile.  The  sheer  profile 
generated  using  the  parameters  displayed  in  Fig.  9  is  plotted  in  Fig.  10. 
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Fig.  9.  The  prompt  box_ for  sheer  profile 
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Fig.  10.  Sheer  profile  curve 


Keel  Rise  Curve 


The  keel  rise  curve  describes  the  shape  of  the  keel  of  the  ship  as  viewed  from  the  side.  To 
generate  this  curve,  HullGen  must  first  calculate  the  section  coefficient  curve,  which  is 
defined  as  the  local  area  of  the  station  divided  by  the  local  beam  on  the  design  waterline 
and  local  draft  as  follows: 


C7(x) 


£4(*)*CX 

DWL(x) 


where  SA(x)  is  local  section  area  curve,  Cx  is  maximum  section  area  coefficient,  and 
DWL(x)  is  local  beam.  When  calculating  this  curve,  the  user  assumes  an  initial  non- 
dimensional  constant  draft  of  value  1 .  The  user  can  then  enter  in  their  desired  position 
where  the  keel  will  begin  to  rise,  as  well  as  a  desired  section  coefficient  at  Station  20, 
Cijo  ,  which  is  shown  in  the  following  prompt  box: 
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Fig.  11.  The  prompt  box  for  keel  rise  curve 

Fig.  12  shows  two  section  coefficient  curves.  The  lower  curve  in  Fig.  12  represents  the 
section  coefficient  curve  if  no  keel  rise  is  allowed.  The  upper  curve  in  Fig.  12  is  faired 
section  coefficient  curve  with  keel  rise.  It  should  be  noted  that  the  upper  curve  is  faired 
into  the  section  coefficient  curve  at  the  station  of  keel  rise  by  matching  the  ordinate  and 
the  slope  of  the  curve  at  this  point.  The  curve  goes  through  the  selected  section 
coefficient  at  Station  20.  These  three  boundary  conditions  are  used  to  generate  a  2nd  order 
polynomial  that  is  plotted  as  an  upper  curve  in  Fig.  12. 


Fig.  12.  Section  coefficient  curve 

With  this  new  section  coefficient  curve,  the  draft  at  each  position  can  be  recalculated  as 
follows: 


m(y\_  SA{x)*Cx 
^  '  DWL(x)*Cj(x) 

which  can  be  used  to  generate  the  keel  rise  curve  shown  in  Fig.  13. 
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Fig.  13.  Keel  rise  curve 


Flat  of  Bottom  Curve 


The  flat  of  bottom  curve  provides  the  user  with  an  outline  of  which  portions  of  the 
bottom  of  the  ship  is  flat  or  completely  horizontal.  It  also  determines  which  of  the  side 
hull  sections  can  be  created  with  a  polynomial  representation  or  a  full  section 
representation  with  a  radius  of  curvature. 

The  curve  is  split  into  five  sections: 

I.  Xj  <  x  <  1.0 

This  section  represents  the  forward-most  section  of  the  flat  of  bottom  curve,  x, , 
the  start  of  the  flat  bottom  curve  and  the  slope  at  xx  are  user  defined.  The  y 
value  forward  of  x,  is  equal  to  the  halfsiding  selected.  The  halfsiding  is  faired 
to  zero  forward  to  station  0. 

II.  X2  <  X  <  Xj 

This  curve  is  represented  by  a  3rd  order  polynomial  defined  by  the  v,  at  x, , 
y2  at  x2 ,  and  the  slopes  at  these  points.  The  slope  at  x2  is  calculated  and  is  not 
user  specified. 

III.  x3  <  x  <  x2 

The  location  of  x2  and  x3  are  user  defined  and  thus  determine  the  extent  over 
which  full  sections  will  be  used  in  the  body  plan.  Thus  y2  and  y3 ,  and  the 
slopes  of  the  curve  at  these  points  are  calculated.  For  the  full  section,  we  have 
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and 


y  =  DWL{x) 


R 

Bx!2 


{DWL{x)~  SAjx)  *  Cx  )11/2 
1-^/4 


where 

Cx  =  maximum  section  coefficient 

Bx  =  maximum  section  beam 

Hx  =  maximum  section  draft 

DWL(x)  =  load  waterline  offset  as  a  function  of  x 

SA(x)  =  section  area  as  a  function  of  x 


IV.  x4  <  X  <  x3 

Similar  to  Section  II  above. 


V.  - 1  <  x  <x4 

Similar  to  Section  I  above. 

The  parameters  used  to  determine  the  flat  bottom  curve  are  defined  in  the  prompt  box 
shown  in  Fig.  14.  The  flat  bottom  curve  can  then  be  obtained  by  combine  five  curves 
described  above.  Fig.  15  shows  the  flat  bottom  curve  and  load  waterline  curve, 
respectively. 

ft  Flat  Bottom  Cuiirve 

Longitudinal  Location  or  the  Start  of  FBC  (XI ) 

O.fl 

Longitudinal  Location  of  the  End  of  Full  Section  (X2) 

0.25 

Longitudinal  Location  of  the  Start  of  Full  Section  (X3) 

-0.425 

Longitudinal  Location  of  the  End  of  FBC  (X4) 

-0.65 

S  lope  Value  at  the  Start  of  FBC  (S 1 ) 

0 

S  lope  Value  at  the  End  of  FBC  (S4) 

0 

Half  siding  CoefT. 

05 

[  OK  j  |  Cancel  j 


Fig.  14.  Prompt  box  for  flat  bottom  curve 
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Fig.  15.  Flat  bottom  curve  and  load  waterline  curve 


Slope  Curves 


In  order  to  assure  reasonable  longitudinal  fairness  in  the  hull  form  being  produced, 
controls  are  placed  on  the  longitudinal  variation  of  the  slopes  that  define  the  section 
shapes  at  deck  at  edge,  design  waterline  and  keel,  thus  three  slope  curves:  keel  slope, 
waterline  slope,  and  deck  slope. 

All  three  curves  are  implemented  in  the  same  manner:  two  3rd  order  polynomials 
connected  with  a  flat  middle  section.  For  each  polynomial,  the  user  inputs  the  necessary 
boundary  conditions  for  the  end  points  of  each  3rd  order  polynomial.  The  user  also  inputs 
the  length  of  the  flat  middle  section,  as  well  as  the  center  point  at  which  the  middle  of  the 
flat  section  occurs.  Each  curve  is  generated  using  a  very  similar  method  to  that  of  the 
deck  at  edge  curve.  Fig.  16  shows  the  prompt  box  for  these  three  slope  curves,  and  Fig. 

1 7  shows  the  three  slope  curves  generated  using  the  parameters  specified  in  the  prompt 
box  shown  in  Fig.  16. 


15 


Center  for  Innovation  in  Ship  Design 


V 


Knowledge  People  Innovation 

Naval  Surface  Warfare  Center  Carderock  Division 

Hydrodynamic  Design  Optimization  Tool 


Section  Slopes  at  Keel 

Section  Slopes  at  DWL 

Section  Slopes  at  Deck 

Slope  Value  at  Station  0 

Slope  Value  at  Station  0 

Slope  Value  at  Station  0 

85.092 

90 

65 

Slope  Rate  at  Station  0 

Slope  Rate  at  Station  0 

Slope  Rate  at  Station  0 

200 

45 

0 

Location  of  Flat  Middle  Section  (Xmid) 

Location  of  Flat  Middle  Section  (Xmid) 

Location  of  Flat  Middle  Section  (Xmid) 

0 

hi.  1023 

lf-0.1025 

Lenqth  of  Flat  Middle  Sections  (Lmid) 

Lenqth  of  Flat  Middle  Sections  (Lmid) 

Lenqth  of  Flat  Middle  Sections  (Lmid) 

0.1 

0.25 

0.1 

Slope  Value  at  Forward  End  of  Flat  Middle  Sections 

Slope  Value  at  Forward  End  of  Flat  Middle  Sections 

90 

Slope  Value  at  Forward  End  of  Flat  Middle  Sections 

4.0 

90 

Slope  Rate  at  Forward  End  of  Flat  Middle  Sections 

Slope  Rate  at  Forward  End  of  Flat  Middle  Sections 

Slope  Rate  at  Forward  End  of  Flat  Middle  Sections 

0.7342 

-12 

o 

Slope  Value  at  Station  20 

Slope  Value  at  Station  20 

Slope  Value  at  Station  20 

20 

48 

80 

Slope  Rate  at  Station  20 

Slope  Rate  at  Station  20 

Slope  Rate  at  Station  20 

0 

45 

0 

(  OK  |  | Cancel) 

[  OK  J  (Cancel) 

|  OK  J  (Cancel) 

Fig.  16.  Prompt  box  for  section  slope  curves 


Fig.  17  Section  slopes  vs  longitudinal  location 


Body  Plan 


The  following  curves  were  designed  with  a  flipped  axis,  i.e.  the  x-axis  is  vertical  and  the 
y-axis  is  horizontal.  This  allows  for  HullGen  to  create  vertical  hull  sides  that  would 
otherwise  be  impossible  to  generate. 

Fine  Sections  (Below  Load  Waterline) 


The  fine  section  is  one  method  of  creating  a  hull  side  underwater.  It  is  narrow  and  has  a 
more  triangular  shape  than  the  full  section,  which  will  be  discussed  later  on.  This  curve 
is  normalized  using  the  local  beam  and  the  local  draft. 
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The  following  boundary  conditions  have  to  be  satisfied: 


(1)  x  =  0.0,  y  =  HS 

(2)  x  =  0.0,  y'  =  Dead  Rise  Slope 

(3)  x  =  l,  y  =  l.O 

(4)  x  =  \,y'  =  slopeDWL 

(5)  \lydx  =  Area 


(Half  siding  (HS)  is  the  normalized  half 
siding  of  the  keel  plate.  This  offset  can 
also  be  the  result  of  a  flat  bottom 
requirement) 

(Dead  rise  curve  is  user  defined  from  the 
keel  slope  curves.  This  value  must  be  great 
than  0) 

(By  definition) 

(Slope  of  the  section  at  the  design  waterline 
as  defined  by  the  slope  curve) 

(Required  area  as  defined  by  the  sectional 
area  curve) 


Instead  of  using  a  regular  4th  order  polynomial  for  above  boundary  conditions,  the  the 
following  polynomial  better  fit  the  required  shape  (Bjorklund  and  Fuller,  1976): 

y  =  AqX  +  Axx  +  A2  +  ^4j(x  +  l)  +  A^(x  +  k )) 

The  constant  k  was  an  arbitrary  value  to  avoid  undefined  results,  which  may  occur  when 
dealing  with  derivatives.  We  used  the  value  k  =  0.0001  in  the  HullGen  2011  program. 
Fig.  18  shows  the  fine  body  section  below  load  waterline  in  a  normalized  coordinate 
system. 


Fig.  18  Fine  body  section  (below  waterline) 
Fine  Sections  (Above  Load  Waterline) 
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The  coordinate  system  is  the  same  as  the  ones  used  for  the  hull  shape.  We  can  also  use 
fine  sections  to  create  the  hull  shape  above  water.  These  are  calculated  using  the  beam 
width,  as  well  as  using  the  difference  between  the  sheer  profile  and  the  draft.  Thus,  this 
curve  is  dimensionalized,  but  it  has  also  been  translated  by  the  draft. 


The  boundary  conditions  for  this  curve  are: 


(1) 

O 

II 

cT 

II 

X 

(Results  of  normalization) 

(2) 

x  =  0,  y'  =  slope DWL 

(This  slope  is  the  same  as  used  at  the  DWL 
for  the  section  below  the  DWL) 

(3) 

x  =  xl,y=yl 

(Local  depth  and  beam  obtained) 

(4) 

x  =  0,  y'  =  slope  Deck  at  edge 

(Flare  slope  is  derived  from  the  slope  curve 
of  the  section  at  the  design  waterline) 

With  these  boundary  conditions,  the  above  water  hull  side  can  be  represented  using  a  3rd 
order  polynomial.  Fig.  19  shows  the  fine  body  section  above  load  waterline  with  actual 
values. 


Fig.  19  Fine  body  section  (above  waterline) 


Full  Sections 


Full  sections  are  characterized  by  having  a  flat  bottom  and  vertical  sidewall,  connected 
by  a  quarter  circle  with  a  constant  radius.  To  create  this  curve,  we  cannot  use  non- 
dimensionalized  coordinates  due  to  distortions  of  the  curve  creating  an  elongated  hull. 

HullGen  first  calculates  the  radius  of  the  hull  using  the  formula: 

R  =  sqrt((H(x)  *  B(x)  -  Area(x))/(  1  —  ir/4)) 

where  H(x)  is  actual  draft  at  particular  station  (feet),  B(x)  actual  offset  at  particular  station 
(feet),  and  Area(x)  actual  area  of  station  as  derived  from  section  area  curve.  With  R 
known,  the  offsets  are  calculated  as: 
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y  =  B(x )  -R  +  V2  *7?*x-x2,  for  0  <  x  <  R 


It  is  just  as  easy  to  create  the  above  water  hull  sides  using  the  same  algorithm  as  these 
with  the  fine  section  hulls;  the  result  will  be  a  flat  line  to  correspond  with  the  flat  side  of 
the  hull.  The  full  body  section  based  on  previous  curves  is  shown  in  Fig.  20. 

40  - , - , - , - , - , - , - 

35  - 

30  - 

25  - 

20  - 

15  - 

1 

10  - 

5  - 

0 - 1 - 1 - 1 - 1 - - ' - 

0  10  20  30  40  50  60  70 


Fig.  20  Full  body  section 


Final  Body  Profile 


With  all  of  these  polynomials  generated,  HullGen  combines  them  to  create  a  body 
profile.  The  user  has  options.  They  can  plot  a  3D  hull  surface  or  body  profiles  by 
combing  the  polynomials  below  and  above  load  waterline. 

To  ensure  that  both  the  methods  developed  by  Bjorklund  and  Fuller  (1976)  were 
accurate,  as  well  as  verify  our  own  code,  we  used  sample  parameters  presented  in  the 
previous  sections  to  create  an  aircraft  carrier  type  ship.  The  polynomial  representation  of 
the  carrier  body  profile  is  shown  in  Fig.  21. 
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Body  Profile 


Fig.  21.  Body  profile  of  a  carrier 

HullGen  has  also  been  used  to  create  a  destroyer/cruiser  type  hull  with  the  data  file 
shown  in  appendix.  The  body  profile  is  shown  in  Fig.  22. 


Body  Profile 


Fig.  22.  Body  profile  of  a  destroyer/cruiser  type  hull 


Conclusions 
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A  rapid  ship  hull  form  generation  tool  is  developed  and  validated.  The  hull  lines/curves 
in  this  tool  are  described  by  polynomials  of  various  orders  and  the  polynomials  are 
determined  using  the  hull  form  parameters.  The  hull  generation  tool  is  implemented  using 
Matlab  to  enable  graphic  user  interface.  The  present  tool  is  able  to  generate  the  same  hull 
forms  as  the  sample  hull  forms  given  by  Bjorklund  &  Fuller  (1976)  with  given 
parameters.  This  rapid  hull  form  generation  tool  describes  hull  form  using  naval 
architect’s  language;  generates  hull  form  from  scratch  in  terms  of  hull  form  parameters; 
and  establishes  the  link  between  the  hull  form  and  form  parameters.  It  can  be  used  to 
generate  the  hull  form  for  hydrostatic  analysis  and  hydrodynamic  design  optimization  in 
the  early  stage  of  ship  design. 


Future  work  will  be  centered  on  comparing  the  seakeeping  results  obtained  using  the 
exact  stations  and  the  approximate  stations  generated  using  HullGen.  The  HullGen  will 
also  be  used  to  generate  various  hull  forms  and  for  optimization. 
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Appendix 


The  following  table  shows  the  sample  input  data  file  for  generating  the  body  profile  of  a 
destroyer/cruiser  type  hull  shown  in  Fig.  22,  and  the  body  profile  of  a  carrier  shown  in 
Fig.  21. 


Input  data  for  a  destroyer/cruiser 
shown  in  Fig.  22 

Input  data  for  a  carrier 
shown  in  Fig.  21 

MainParameters 

MainParameters 

numStations 

numStations 

20 

20 

Length 

Length 

600 

925 

Beam 

Beam 

70 

127 

Draft 

Draft 

23 

35.5 

maxCoef f 

maxCoef f 

0.82 

0.98 

numSegments 

numSegments 

20 

40 

SectionAreaCurve 

SectionAreaCurve 

A0 

A0 

0.0 

0.001 

SO 

SO 

-1 . 123 

LO 

o 

i 

Xmax 

Xmax 

-0 . 052 

-0 . 1 

Smax 

Smax 

0 

0 

Lmid 

Lmid 

0 

0.001 

A2  0 

A2  0 

0.04 

0 . 05 

S20 

S20 

0.9718 

1 . 0 
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Cp 

Cp 

0.62 

0.63 

LCB 

LCB 

-1 

-3 

DesignWaterline 

DesignWaterline 

Y0 

Y0 

0.005 

0.005 

SO 

SO 

i 

M 

o 

00 

-1 . 

Xmax 

Xmax 

-0.052 

-0 . 11 

Smax 

Smax 

-0.025 

0 

Lmid 

Lmid 

0 

0 . 07 

Y20 

Y20 

0.55 

0.4063 

S20 

S20 

1 

1.65 

Cwp 

Cwp 

0.76 

0.7514 

LCF 

LCF 

LO 

LO 

i 

-6.0 

DeckEdge 

DeckEdge 

Y0 

Y0 

0.436 

0.75 

SO 

SO 

-1.0649 

LO 

CM 

O 

1 

Xmid 

Xmid 

-0 . 1025 

-0 .1025 

Lmid 

Lmid 

0.22 

0.22 

Y10 

Y10 

0.989 

1.0 

S10 

S10 

-0.06103 

0 
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Y20 

Y20 

0.6592 

0.9 

S20 

S20 

0 . 9427 

0.35 

Sheer 

Sheer 

DO 

DO 

60 

72 

DIO 

D10 

43 

72 

D20 

D20 

44 

72 

Drafts 

Drafts 

keelPos 

keelPos 

-0.15 

LO 

LO 

o 

i 

coef f 20 

coef f 20 

0.65 

0.65 

BottomSlope 

BottomSlope 

TO 

TO 

70 

85.092 

SO 

SO 

240 

200 

Xmid 

Xmid 

-0 . 052 

0 

Lmid 

Lmid 

0 

0 . 1 

T10 

T10 

4 

4.0 

S10 

S10 

6.547 

0.7342 

T20 

T20 

1 . 4474 

20 

S20 

S20 

20 

0 

water line SI ope 

water lineS lope 

TO 

TO 
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88 

90 

SO 

SO 

35 

45 

Xmid 

Xmid 

-0 . 052 

-0.1023 

Lmid 

Lmid 

0 

0.25 

T10 

T10 

90 

90 

S10 

S10 

0 

-12 

T20 

T20 

77 . 968 

48 

S20 

S20 

20 

45 

deckSlope 

deckSlope 

TO 

TO 

69 

65 

SO 

SO 

-10 

0 

Xmid 

Xmid 

-0.052 

-0.1025 

Lmid 

Lmid 

0 

0 . 1 

T10 

T10 

89 

90 

S10 

S10 

-10 

0 

T20 

T20 

84 

80 

S20 

S20 

0 

0 

f latBottom 

f latBottom 

XI 

XI 

LO 
\ — 1 

o 

1 

0.8 

X2 

X2 
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LO 
\ — 1 

o 

1 

0.25 

X3 

X3 

-0.15 

-0.425 

X4 

X4 

-0 . 15 

-0.65 

SI 

SI 

0 

0 

S4 

S4 

0 

0 

HS 

HS 

0.5 

0.5 
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